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ABSTRACT 

Number Counts of galaxies in two GALEX bands (FUV: 1530A and NUV: 
2310A, both in AB magnitudes) are reported. They provide for the first time in 
the literature homogeneously calibrated number counts of UV galaxies covering 
continuously a very wide range of UV magnitude (14 - 23.8). Both the FUV 
and NUV counts are inconsistent with a non-evolution model, while they are in 
good agreement with evolution models (essentially luminosity evolution) derived 
from the high-z UV luminosity functions of Arnouts et al. (2004). It is found 
that the contribution from galaxies detected by GALEX to the UV background 
is 0.68±0.10 nW m'^ sr'^ at 1530A and 0.99±0.15 nW m'^ sr"^ at 2310A. 
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These are 66±9% and 44±6% of the total contributions of galaxies to the the 
UV background at 1530A (1.03±0.15 nW m'^ sr'^) and at 2310A (2.25±0.32 
nW sr^^), respectively, as estimated using the evolution models. Galaxy 
counts and star counts in 7 regions, each contains a few dcg^ GALEX coverage 
in an area of up to ~ 30 deg^, arc compared with each other to study the region- 
by-region variance. This shows that for the galaxy counts the cosmic variance is 
comparable to the net error due to other uncertainties. The star counts increase 
with decreasing absolute Galactic latitude |b|. 

Subject headings: ultraviolet: galaxies - galaxies: photometry - galaxies: active 
- galaxies: evolution - galaxies: stellar content - stars: formation 

1. Introduction 

The number counts of UV galaxies as a function of the UV magnitude provide direct 
measures of the density and the evolution of star forming galaxies. They also give strong 
constraints on the brightness of the cosmic UV background. In the literature, the UV 
number counts have only been measured in a few small sky areas and in some discontinued 
UV magnitude ranges by several inhomogeneously calibrated UV surveys (Milliard et al. 
1992; Deharveng et al. 1994; Gardner et al. 2000; Sasseen et al. 2002; Iglesias-Paramo et 
al. 2004). The Galaxy Evolution Explorer (GALEX) has opened a new era of extra-galactic 
UV astronomy. It is surveying the UV sky with unprecedented efficiency and very high 
sensitivity. With a pyramid-like survey structure (Martin et al 2004), GALEX presents a 
uniform and complete picture of galaxy evolution in the UV domain since z~ 1. In this 
paper, we carried out galaxy number counts in the two GALEX bands: the FUV (1530A) 
and the NUV (2310A). We have the following science goals: (1) constraining the evolution 
models of UV galaxies; (2) evaluating the cosmic variance of UV galaxies; (3) estimating 
the contribution from galaxies to the cosmic background radiation in the two UV bands; (4) 
providing the calibrations for the UV luminosity functions. 

2. Data 

Three data sets taken from the GALEX database were analyzed in this work, yielding 
17174 FUV galaxies and 41512 NUV galaxies in the final samples for the counts (Table 1). 
All the magnitudes are in the AB system. The calibrations (1R0.2 calibration) have errors of 
~ 10% for both FUV and NUV (Morrissey et al. 2004). All of the magnitudes are corrected 
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for the Galactic extinction using the Schlegel et al. (1998) reddening map and the Galactic 
extinction curve of CardeUi et al. (1989). 

The main data set includes 36 Medium-depth Survey (MIS) fields taken from the 
GALEX internal data release IR0.2. The exposure time of these fields is in the range of 
1000 - 1700 seconds with the median of ~ 1500 seconds. They all have full coverage in 
the Sloan Digitized Sky Survey (SDSS hereafter) DRl database (Abazajian et al. 2003). 
In order to minimize the artifacts which concentrate near the periphery of the field of view 
(Morrissey et al. 2004), we include only sources within 0.45 deg radius from the field center, 
corresponding to a sky coverage of 0.64 deg^ for each field. The source lists are products 
of the GALEX pipeline, which uses a combination of the programs SExtractor (Bertin & 
Arnouts 1996) and Poissonbg, a new program written for GALEX data, to detect and obtain 
photometry for sources in GALEX images. Poissonbg determines a background map for each 
GALEX image with a method very similar to that used in SExtractor itself except that the 
background is estimated using the mean and counts due to sources are iteratively clipped 
out using the full Poisson distribution, rather than assuming Gaussian statistics which is not 
appropriate for low counts rate found commonly in GALEX images. Poissonbg also creates 
a detection threshold map for a user-specified probability using the background image and 
the full Poisson distribution. The original image is divided by the detection threshold map 
and given as the detection image to SExtractor while the SExtractor measurements are done 
on the background-subtracted image. The GALEX magnitude is taken from the MAG auto 
of the SExtractor. Key SExtractor/Poissonbg parameters used by the Pipeline are hsted 
in Table Al (electronic table). The second data set contains 3 Deep Survey (DIS) fields 
{Texp =25536s, 23591s, 5607s), all in the Spitzer First Look Survey (FLS) area. These fields 
are also covered by SDSS DRl. In these DIS fields the source confusion becomes significant. 
Hence fine tunings in the SExtractor/Poissonbg parameters have been applied (see electronic 
Table A2). A slight dimming (on the order of 0.1 - 0.2 mag) of the fluxes at magnitudes 
fainter than 23 was found in the later stage of data analysis, resulting from a small (~ 2.5%) 
positive bias of our present sky background estimator. This was corrected empirically ac- 
cording to results of the artificial source simulations (see Section 3) . The minimum distance 
among the 3 DIS fields is 0.86 deg. In order to avoid source duphcation, sources were in- 
cluded only when they are within 0.43 deg radius from the field center, corresponding to a 
coverage of 0.58 deg^ per DIS field. The third data set consists of 95 bright GALEX galaxies 
(NUV < 16) selected by Buat et al. (2004) in the study of the FIR emission of UV galaxies. 
The characterization of these galaxies can be found in Buat et al. (2004). 

Matches with the SDSS catalogs were carried out for the MIS and DIS sources with a 
search radius of 4". Less than 3% of FUV and NUV sources brighter than 22.5 mag (extinc- 
tion corrected) have no SDSS counterparts. These are predominantly spurious sources due to 
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artifacts such as pieces of shredded bright extended sources. However, optical counterparts 
of r > 22 of fainter GALEX sources may be missed by SDSS. In the magnitude range of 22.5 
and 23.8, about 10% and 20% of FUV and NUV sources do not have SDSS counterparts, 
respectively. Detailed inspections of the GALEX images showed that most of these sources 
are real. The SDSS star/galaxy classifications were adopted for MIS/DIS sources with SDSS 
counterparts. Sources brighter than 22.5 mag and without SDSS counterparts are deemed 
spurious and dropped off from the sample. Sources with UV magnitudes in the range of 22.5 
and 23.8 mag and without SDSS counterparts were included as galaxies. 



3. Bias Corrections 

Incompleteness: We adopt the Monte Carlo simulation algorithm developed by Small 
et al. (1995) in assessing the incompleteness of GALEX catalogs. Artificial sources of a 
given UV magnitude were added randomly to a GALEX image. Then the same source 
extractor which produced the original catalog from the same image was applied and the 
results were checked for non-detections of the added sources. This gives an estimate of the 
incompleteness. In order to reproduce the real point spread function (PSF) and the photon 
noise, the artificial sources were created by dimming the bright sources extracted from the 
same image, and the photon counts in the region where an artificial source is added were 
randomized according to the Poisson probability function. The counts in each field were 
truncated at the magnitude where the incompleteness is becoming larger than 20%. In 
magnitude bins brighter than the completeness cut-off, the error due to the incompleteness 
correction is estimated conservatively to be 50% of the correction. 

Spurious Sources: The spurious source fraction as function of the non-dereddened apparent 
magnitude was estimated by counting single-detection sources in regions with multipule 
coverages ('repeatability' method). The effect of incompleteness was accounted for in the 
analysis. The error is estimated to be 50% of the correction. 

Star/Galaxy Classification Errors and Contamination of AGNs: The biases due to 
mis-classifications of sources by SDSS and due to AGN contamination are accessed using 
the results of photo-z processing of GALEX sources in the MIS fields. By fitting fiuxes in 7 
bands (5 SDSS bands plus 2 GALEX bands) with templates of stars, galaxies and AGNs, the 
photo-z processing makes the best use of the information in the SEDs. Although there are 
still substantial uncertainties for individual sources, the statistics derived from these results 
are robust. Galaxy counts in a given magnitude bin were derived using the formula 

Cgal — Cext X Qext + Cpnt X Qpnt = Cgxt ^ 9 gal (1) 
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Fig. 1. — Left: FUV (1530A) number counts of galaxies. Right: NUV (2310A) number 
counts of galaxies. Number counts in other UV bands are converted to counts in the GALEX 
bands according to an assumed UV slope of (3 = —0.8. 



Table 1. FUV and NUV Number Counts of Galaxies 
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0.016 ± 0.012 


1.00 


1.00 


4 


615.00 



Note. — Table 1 is published in its entirety in the electronic edition of The Astrophysical Journal 
Letters. A portion is shown here for guidance regarding its form and content. 

^/ = (1 - /sp«r)/(l - fincomp), whcrc fspur IS the corrcctiou for the spurious sources, and fincomp the 
correction for the incompleteness. 
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where Cgxt is the number counts of SDSS extended sources, Qext the fraction of galaxies among 
SDSS extended sources, Cpnt the number counts of SDSS point sources, gpnt the fraction of 
galaxies among these sources, and ggai = gext + Qpnt x Cpnt/cext- Both g^xt and gpnt were taken 
from the photo-z results. A 5% error was assigned to the MIS and DIS counts in all FUV 
and NUV magnitude bins to take into account the uncertainty due to mis-classifications. 

Errors in Foreground Extinction Correction: The following errors can occur to the 

extinction correction: (1) incorrect UV extinction law in average; (2) region-to-region varia- 
tions in the UV extinction law; (3) insufficient angular resolution of the Schlegel map. Based 
on results of an extensive test involving large samples of low and high extinction regions, a 
conservative 10 percent extinction correction error was assigned to the counts. 

Eddington Bias: The Eddington bias (Eddington 1913) is corrected using the algorithm 
developed by Hogg & Turner (1998, Eq.(4)), assuming the slope of the counts near the 
detection limit to be p = 1. 

4. Results 

The number counts of galaxies in the FUV and NUV bands (Table 1) cover continuously 
the magnitude range of 14.0 - 23.8 for FUV and 14.0 - 23.6 for NUV. The error bars include 
not only Poisson errors, errors due to incompleteness and spurious sources corrections, and 
errors due to source classifications and extinction correction, but also errors due to cosmic 
variance calculated using the formulism of Glazebrook et al. (1994), which are in general on 
the same order of Poisson errors. Note that the cosmic variance might be underestimated 
by the assumption that individual GALEX fields are independent with each other. In Fig.l 
the results are compared with UV counts taken from the literature, and with predictions of 
evolution models. Both the FUV counts and the NUV counts are slightly lower than the 
FAUST and FOCA counts (Dcharveng et al. 1994; Milhard ct al. 1992; Iglesias-Paramo et 
al. 2004) in bins brighter than 21 mag. This is likely to be due to a calibration difi^crcnce 
between GALEX and FAUST/FOG A (see also Wyder et al. 2004). It should also be noted 
that counts of bright UV galaxies by Deharveng et al. (1994) and Iglesias-Paramo et al. 
(2004) are biased to higher values by galaxy clusters. The XMM- Newton 0M-UVW2 counts 
(Sasseen et al. 2002) appear to be incomplete in bins fainter than 21 mag, where they 
decrease with increasing magnitude. The GALEX counts are consistent with the HST counts 
(Gardner et al. 2000) in the magnitude bins where the two overlap. 

Both the FUV and NUV counts are inconsistent with predictions by a non-evolution 
model, which was constructed using the local luminosity functions in the FUV and NUV 



- 7- 




17 18 19 20 21 22 17 18 19 20 21 22 
FUV NUV 

Fig. 2. — FUV and NUV number counts of galaxies and stars in 7 regions with MIS and 
DIS coverage. 
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(Wyder et al. 2004). The K-correction was estimated using an SED based on the SB4 
class spectrum of Kinney et al. (1996) at A > 1200A, a dropping spectrum (by a factor of 
~ 2) between 1200 — lOOOA, and a sharp cut-off shortward of 912A. The SB4 spectrum has 
f3 = —0.8. This is close to the value calculated from the ratio of the FUV and NUV luminosity 
densities in the local universe by Wyder et al. (2004), after allowing for the difference in the 
calibration zero points (0.03 mag for FUV and -0.10 mag for NUV, Morrissey et al. 2004) 
used in these two works. The high redshift UV luminosity functions of Arnouts et al. (2004), 
which indicate a UV luminosity density evolution of (1 -|- 2;)^-^='=° ''' to z ~ 1 (Schiminovich et 
al. 2004), can be best fitted by an analytic model with oc (1 -|- z)"^-^, 0* oc (1 -|- z)~^-^ and 
a oc (1 -|- z)~^^^. In the expression for the integrated luminosity density of different redshift, 
Pl(z) = 0*(z)L^(z)r(2 — a(z)), the evolution of 0^, and that of a cancel with each other for 
< a < 2, leading to Pl(z) ~ L*(z). Therefore this is essentially a luminosity evolution 
model. Predictions by two models, both assuming the same evolution parameters as given 
above but using different UV SEDs, arc plotted in Fig.l: Model 1 assumes the same SED 
as that of the no evolution model, and the SED assumed in Model II is otherwise the same 
except for a flat spectrum between 1200 — lOOOA. These two models should embrace the 
uncertainties in the SED of wavelengths shorter than the Lya due both to the Lya absorption 
and dust extinction (e.g. Buat et al. 2002), which affects signiflcantly the K-corrections in 
GALEX bands. In the bins covered by the GALEX counts, the two evolution models give 
almost identical results, both fit the data very well. In fainter bins, the two models show 
significant difference in the FUV band due to different K-corrections. Both models fit the 
HST NUV counts well, given the substantial uncertainties of the data, but are slightly lower 
than the HST FUV counts. 

Using the mean of the integral counts predicted by Model I and Model II, and the beam 
size of the DIS FUV maps (FWHM=5.4") and that of DIS NUV maps (5.6"), we estimate 
that the 40-beam-per-source confusion limits of the FUV and NUV maps are FUV=25.3 and 
NUV=24.0, respectively. The major contribution to the incompleteness at the boundary of 
the NUV counts (NUV=23.6) is from the source confusion, indicating the NUV counts being 
confusion limited. The confusion hmit of the DIS FUV maps, on the other hand, is signifl- 
cantly deeper than the last bin of the FUV counts. The contribution from galaxies detected 
by GALEX to the UV background is 0.68±0.10 nW m'^ sr-^ at 1530A and 0.99±0.15 nW 
m~^ sr~^ at 2310A. Extrapolating the counts using the above two evolution models and inte- 
grating down to flux=0, the total extragalactic UV background due to galaxies is 1.03±0.15 
nW m~^ sr~^ at 1530A and 2.25±0.32 nW m ^ sr at 2310A, respectively. The errors 
include both model uncertainties and GALEX calibration uncertainties. Our results are sig- 
niflcantly lower than the results of Gardner et al. (2000) who found that the contributions 
from resolved sources to the UV background are 3.9loi m~^ sr' -1 and S.Gtoi nW m 
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sr~^ at I595A and 2365A, respectively. Our results are in good agreement with that of 
Armand et al. (1994), who found that the contribution from galaxies to the UV background 
at 2OOOA is 40 - 140 photons cm~^ sr~^, corresponding to 0.80 - 2.78 nW m~^ sr~^. 

In Fig. 2 FUV and NUV number counts of galaxies and stars in 7 regions with MIS and 
DIS coverage are presented in 4 panels. For galaxy counts the error bars do not include the 
cosmic variance in these plots. It appears that differences between galaxy counts in different 
regions are consistent with the error bars, indicating that the cosmic variance is not very 
significant for UV surveys covering more than a few square degrees. The counts of stars 
show a significant trend in the sense that the less the absolute Galactic latitude (|b|), the 
higher the counts. The NUV counts of stars are rather flat, higher than galaxy counts in 
bins brighter than ~ 21 mag. The FUV counts of stars are lower than the counts of galaxies 
in bins fainter than ~ 18 mag. These data are compared with predictions of a star counts 
model. The model closely follows the Bahcall-Soneira model (Bahcall 1986) with an updated 
parameterization of the Galactic spheroid from Gould et al. (1998). The stellar luminosity 
function of Wielen et al. (1983) is used. An additional disk white dwarfs population (scale 
height of 275 kpc) is included, with a luminosity function drawn from Liebert et al. (1988). 
The model generally reproduces the trends of the observed star counts, though the predicted 
counts in both bands are slightly flatter than the observations. 
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